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ABSTRACT 

Our 50 ks Chandra observation of the small radio supernova remnant (SNR) Gl .9+0.3 shows a complete shell 
structure with strong bilateral symmetry, about 100" in diameter. The radio morphology is also shell-like, but 
only about 84" in diameter, based on observations made in 1985. We attribute the size difference to expansion 
between 1985 and our Chandra observations of 2007. Expansion is confirmed in comparing radio images from 
1985 and 2008. We deduce that Gl. 9+0.3 is of order 100 years old - the youngest supernova remnant in 
the Galaxy. Based on a very high absorbing column density of 5.5 x 10 22 cm -2 , we place Gl. 9+0.3 near the 
Galactic Center, at a distance of about 8.5 kpc, where the mean remnant radius would be about 2 pc, and the 
required expansion speed about 14,000 km s" 1 . The X-ray spectrum is featureless and well-described by the 
exponentially cut off synchrotron model srcut . With the radio flux at 1 GHz fixed at 0.9 Jy, we find a spectral 
index of 0.65 and a rolloff frequency of 1.4 x 10 18 Hz. The implied characteristic rolloff electron energy of 
about 94(5/10 /iG) -1 / 2 TeV is the highest ever reported for a shell supernova remnant. It can easily be reached 
by standard diffusive shock acceleration, given the very high shock velocities; it can be well described by 
either age-limited or synchrotron-loss-limited acceleration. Not only is Gl. 9+0.3 the youngest known Galactic 
remnant, it is also only the fourth Galactic X-ray synchrotron-dominated shell supernova remnant. 
Subject headings: supernova remnants, X-rays : general — supernova remnants: individual (Gl. 9+0.3) — 
X-rays: ISM 



1. INTRODUCTION 

Estimates for the supernova rat e in the Milky Way Galaxy 
give about 3 per century (e.g., Ivan den Bergh & Tammanl 
Il99l|) . so that we expect 60 supernova remnants (SNRs) 
younger than 2000 yr, while fewer than 10 are confirmed. 
This well-known deficit motivated our program of Chandra 
observations of compact radio remnants, in search of evi- 
dence of youth. Here we report a confirmation of a young 
SNR, but one of remarkable properties: comparison of radio 
and X-ray images separated by 22 years demonstrates signifi- 
cant expansion, allowing the inference of an age of order 100 
years. In addition, its X-ray emission appears to be purely 
synchrotron radiation, making Gl. 9+0.3 the fourth clear-cut 
X-ray- synchrotron- dominated Gal actic SNR, after SN 1006 
dKovamaetalJI 1995k G 347.3-0.5 dKovama et al.ll 19971) . and 
G266.2-1.2 dSlane et alJl200lh . 

Gl. 9+0.3 was identified bv lGreen & Guil (11984 as a po- 
tential young SNR, from a radio image showing a 1.'2- 
diameter shell. This angular size made Gl. 9+0.3 the small- 
est Galactic SNR, with a mean radio surface brightness about 
half that of Tycho's SNR and about 20 times that of SN 1006, 
suggesting an age of order 10 3 yr. Green & Gull (1984) 
point out the general resemblance of Gl. 9+0.3 to Kepler's 
SNR, and mention that it would be below 10 pc in diame- 
ter at any distance in the Galaxy. VLA 4 observations from 
1985 (Green 2004, reproduced in Fig.Q]) show a strong asym- 
metry in the shell at 21 cm, perhaps indicative of an exter- 
nal density gradient. Radio flux measurements with differ- 
ent instruments and frequencies are inconsistent; from VLA 
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FIG. 1.— 1985 radio image of Gl.9+0.3 (Green 2004): VLA at 1.5 GHz. 
Resolution 8.4" x 3.6", position angle 3.5° E of N. Intensities in Jy/beam. 

measurements, iNord et ail (12004]) report a radio spectral in- 
dex a ~ 0.7 (S u oc v~ a ) between 74 and 327 MHz (consistent 
with flux estimates by Gray 1994) and giving a 1 GHz flux 
of 0.9 Jy. This steep a spectrum is typical of historical SNRs 
but rare otherwise. No counterparts are obvious in optical 
or infrared images. The remnant was detected with ASCA 
( Sak ano et al.l l2002|) . but contamination from nearby bright 
sources precluded detailed analysis. 

2. X-RAY IMAGES AND SPECTRA 

We observed Gl.9+0.3 with Chandra using the ACIS-S 
CCD camera (S3 chip) on 2007 February 10 for 24 ks, 
and March 3 for 26 ks. We checked aspect correction 
and created new level- 1 event files appropriate for VFAINT 
mode. No flares occurred during the observation. CTI 
correction was applied and calibration was performed using 
CALDB version 3.4.0. Finally, the datasets were merged and 
weighted response files created. We extracted spectra using 
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Fig. 2. — Chandra image of Gl. 9+0.3, platelet smoothed, covering the 
same area as Figure 1. Colors are counts/ACIS pixel (0.49 // )between 1.5 and 
6 keV. 

thespecextract script, using a background region on the 
S3 chip. About 8,000 counts were obtained from the remnant. 

Figure [2] shows the image over the full energy range, 
platelet- smoothed (method described in Willett 2007). The 
bright ring is nearly circular, though fainter extensions 
("ears") protrude symmetrically on the E and W sides. The 
azimuthal brightness variations are quite different in radio 
and X-rays; the former has a single maximum to the north, 
while the latter shows remarkable bilateral symmetry reminis- 
cent of SN 1006. However, the difference in radio and X-ray 
morphologies is quite unlike SN 1006. The angular extent is 
considerably larger than suggested by the radio image: about 
120" x 90". The bright ring is slightly smaller E-W, about 
100" x 90". No central source is apparent. 

Instead of subtracting background (rendering the statistics 
non-Poisson), we modeled it with a combination of Gaussians 
and power-laws, exponentially cut off on both ends. The spec- 
trum is shown in Figure along with the background model. 
No lines at all are evident. Very high absorption is also appar- 
ent. We fit the data with the srcut model for synchrotron 
emission from a power-law distribution of electrons with a 
high-energy exponential cutoff (Reynolds & Keohane 1999). 
The model has four parameters: absorbing column density 
Nh, 1-GHz radio flux, power-law spectral index a, and rolloff 
frequency z/ ro ii; we fixed the radio flux to 0.9 Jy (though 
brightening may have occurred; Green et al. 2008). Parameter 
estimation was done with Markov Chain Mont e Carlo simula - 
tions, as implemented in the XSPEC package ( Arnaud 1996). 
Constant Bayesian priors were assumed for the three undeter- 
mined parameters. We found Nh = (5.5 ± 0.3) x 10 22 cm -2 , 
a = 0.65 ± 0.02 and i/ roU = 1.4^ x 10 18 Hz (Ai/ roU = 5.8^ 
keV) - one of the highest values ever reported for a SNR. 
(Errors are 90% confidence intervals.) Now the high column 
means that dust scattering will be important, preferentially 
for lower-energy X-rays, hardening the spectrum. Including 
scattering will reduce the inferred column density and rolloff 
frequency somewhat. Detailed modeling will be required for 
quantitative assessment. In any case, evidently we have dis- 
covered the fourth example of a Galactic synchrotron-X-ray 
dominated SNR. 

The striking size difference led us to compare the radio and 
X-ray images in more detail. The bright X-ray ring shows 
indentations and bright features that correlate well with the 
outer extent of the radio image - though at considerably larger 
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• Integrated spectrum of Gl. 9+0.3, binned to a minimum of 25 



counts per channel for display purposes only. Model: srcut (see text). 



Fig. 4. — Comparison of radio images (red) and X-ray image smoothed 
to 10" x A" . Inset: Original 1985 radio image. Main figure: Radio image 
expanded by 16%. Lines: orientation of profile in Fig. 5. 

radius, as can be seen in FigureH] We smoothed the raw X-ray 
image to the resolution of the radio image, and examined an 
E-W profile (Figure 5) chosen to pass through prominent ra- 
dio and X-ray features. Finite resolution should not affect the 
radio/X-ray size ratios. The distance between X-ray peaks, 
about 95", should be a lower limit to the X-ray diameter; if 
the X-rays abruptly turned on at the shock, those locations 
would appear as the outer inflection points, separated by about 
100", a clear upper limit. The outer radio edge was estimated 
from inflection points in the radio profile, not from the peaks, 
as we expect an intrinsically much thinner X-ray than radio 
ring. The radio diameter judged in this way was 84". These 
values give fractional expansions of between 13% and 19% 
of the 1985 size. We adopt an estimate of (16 ± 3)% for the 
expansion in the 22 years since 1985, giving an age estimate 
of about 140 ± 30 yr, or less if deceleration has occurred as is 
likely. We shall argue below that Gl.9+0.3 is near the center 
of the Galaxy; at 8.5 kpc, the shock speed is 14,000 km s -1 . 
While such speeds are seen in optical spectra of supernovae, 
they are unprecedented in supernova remnants. 

We recently obtained VLA Exploratory Time to examine 
Gl.9+0.3 at the present epoch (Green et al. 2008). The VLA 
is currently in C configuration, so we observed at 6 cm to 
achieve comparable resolutions (10" x 4"). Averaged over all 
azimuths, the shell profile peaks at a radius of about 30" in 
1985 and 34.5" in 2008. There are also some morphological 
differences in the images; however, a careful analysis of the 
averaged profiles gives an expansion of (15 ± 2)% (Green et 
al. 2008), consistent with our determination. 

3. DISCUSSION 
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Fig. 5. — Profiles of radio and X-ray images. Outer (red) line: X-ray. 
Dashed green line: Original radio image. Short-dashed blue line: Radio ex- 
panded by 16% about the center of the ring, to bring outer inflection points 
to X-ray maxima. X-ray extensions ("ears") are apparent at about 30% of 
maximum. 

The extremely high absorption toward Gl.9+0.3 makes it 
unlikely that it is much closer than the Galactic Center (GC) 
re gion, at about 8 . 5 kpc. In fact, the Galactic extinction model 
of Marshall et al. ( 2006), based on near infrared 2MASS ex- 
tinction colors, gives an infrared absorption of A Ks = 1.8 m at 
the GC distance and / = 1.87°, b = 0.32°. Assuming the same 
amount of heavy-e lement depletion onto d ust grains as in the 
dust model of Weingart ner & Dra ine ( 2001), about half of the 
absorption can be accounted for by material located in front 
of the GC. However, the very high expansion we infer ar- 
gues against much larger distances than to the Galactic Cen- 
ter, suggesting that the additional absorption is local to the 
source. We assume a nominal distance of 8.5 kpc, at which 
\" = 0.041 pc. Then the mean radius of the bright X-ray ring 
is about 2 pc, with the E and W "ears" at about 2.2 pc. The 
very high expansion velocities would be mitigated at a sub- 
stantially smaller distance, but the extreme absorption and ab- 
sence of any IR counterpart (for instance, from a more local 
dense cloud of some kind) argue strongly against this. 

Our inferred age of about 100 yr makes Gl.9+0.3 the 
youngest known remnant in the Galaxy. The very high ab- 
sorption would have rendered the supernova undetectably 
faint. The average optical absorption Ay foreground to the GC 
exceeds 15 m , and the dust local to the GC region contributes a 
comparable amount, resulting in Ay > 30 m . Even a bright SN 
la with M v ~ -\9 ,n would be fainter than my ~ 25 m . 

Three arguments slightly favor a Type la origin for 
Gl.9+0.3. First, the bilateral symmetry of the X-ray syn- 
chrotron emission suggests interaction with a roughly uniform 
magnetic field, with electron acceleration dependent on the 
obliquity angle between shock normal and B (e.g., Fulbright 
& Reynolds 1990). A stellar- wind magnetic field, as expected 
for a core-collapse (CC) progenitor, should be a tightly wound 
Archimedean spiral, i.e., almost perpendicular to the shock 
normal at all azimuths, removing that possible explanation for 
the nonthermal morphology. However, Tycho, a prototype la 
remnant, does not show bilateral symmetry, nor does the ra- 
dio emission in Gl.9+0.3. Second, the very high expansion 
speeds are slightly more probable for a la supernova. Third, 
there is no indication of a pulsar-wind nebula at the center 
(but a compact object with a soft spectrum like that in Cas A 
could not be seen with this high absorption). The question is 
still open; the position of Gl.9+0.3 near the GC is consistent 
with either type, as both the bulge and disk populations are 
substantial there. 

The high velocities suggest that little deceleration has oc- 
curred. The mass of gas in a sphere of radius 2.2 pc is only 
1.5no Mq for a mean number density no and cosmic abun- 
dances. A model Type la explosion with an exponential ejecta 



profile ( Dwarkada s & Chevalied fT998h achieved the current 
size and shock velocity for standard explosion parameters, 
an age of 100 yr, and no ~ 4 x 10~ 2 cm -3 , characteristic of 
the ISM hot phase. For this density, the extremely low X- 
ray emission measure, combined with the ionization timescale 
r = not of only ~ 10 8 cm -3 s, means that no significant X-ray 
line emission would be expected from shocked ISM. The age 
is so small that very little of the ejecta have been shocked, so 
a contribution to the spectrum from shocked ejecta should be 
small as well. 

Densities might be larger than average in the northern radio- 
bright limb. We partitioned the shell into two sections, the 
northern limb and the remainder, and made joint fits with 
a nonthermal srcut model in both limb sections, adding 
a vp shock model with cosmic abundances in the northern 
limb alone. Thermal emission with emission measure as high 
as 1 .6 M cm -3 for gas with temperature 4 keV and shock age 
2 x 10 9 cm -3 s can be accommodated by the current data; the 
very low ionization age suppresses most X-ray lines. (How- 
ever, the absence of Fe Ka emission confirms that the contin- 
uum we see is nonthermal.) 

The rolloff photon energy is related to the e-folding en- 
ergy of the exponential cutoff in the electron spectrum, E m , 

by E m = 39(hiy ro n) l/2 B~ l /2 TeV ([Reynolds & KeohandfT999l 
with a numerical error corrected), where hv TO \\ is in keV, and 
2?io = B/10 fiG is the post-shock magnetic field. From the 
fitted value hv vo \\ = 5.9 keV, we infer E m = 94 TeV - the high- 
est value ever reported for an integrated spectrum. However, 
with our very high shock velocities, the predictions of stan- 
dard shock acceleration theory can easily accommodate such 
energies. If radiative losses cut off the electron spectrum, 
E m ex Br 1 / 2 , so hv Y0 \\ = 25. lE^B keV is independent of B. 
Standard results (e.g., Reynolds 1998) for a shock with speed 
u% = w/1000 km s -1 give hv m \\ ~ 0.8wg keV, for Bohm diffu- 
sion and a perpendicular shock, and a compression ratio of 4 
(the estimate depends only weakly on these assumptions). We 
estimate u% = 14, so we can easily reach the required value, 
even with considerably less efficient acceleration. An esti- 
mate of the limitation due to finite age can be made assuming 
no deceleration; the expression for acceleration time to energy 
E given in Reynolds (1998) gives r(94 TeV) = 300B^uf yr, 
so that the shock could produce our observed E m in only a few 
years. Likely higher values of E m at the bright limbs would 
then not pose any problem. 

The contrasting X-ray and radio azimuthal brightness pro- 
files imply substantial spectral differences between the north- 
ern and southern hemispheres. Preliminary spectral fits do 
suggest differences, but unexpectedly mainly in a rather than 
z/ ro ii. Efficient shock acceleration is expected to produce 
concave-up curvature in the particle spectrum (hardening to 
higher energies), so that a in srcut is a mean radio-to- 
X-ray spectral index. If deeper X-ray observations confirm 
this form of spectral difference, the result may have signifi- 
cant implications for the understanding of shock acceleration. 
Other radio-X-ray differences, such as the faint "ears" of X- 
ray emission beyond the bright ring on the E and W sides, 
could in principle be due to pre-shock diffusion, or simply re- 
flect lower sensitivity to faint radio emission. Estimates of the 
diffusion length, assuming a mean free path of one gyroradius 
of a 94 TeV electron in a magnetic field of a few /iGauss, are 
roughly consistent with the observed length scale of 5" - 10", 
but more extensive radio and X-ray observations are needed 
to test this explanation. 
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A supernova this recent should produce a significant 
signal in the 1.156 MeV gamma-ray decay line as ra- 
dioactive 44 Ti decays first to 44 Sc and then to 44 Ca, as 
has been reported for Cas A (llyudin et al.l 11994)) . Es- 
timates of the yield of 44 Ti in SNe la range from 8 x 
1(T 6 Mo to 5 x 1(T 5 M© (llwamoto et al.1 Il999h . while 
CC predictions range between 2 x 10~ 5 a nd 2 x 10~ 4 M 
(Thielem ann. Nomoto. & Hashi moto! 1 19961) . For a half-life 
Um = 60.3 ± 1.3 vr (iGorres et al.lll998l) . we predict a flux of 
2.3 x 10- 6 (M(Ti)/l(T 5 Af ) ph cm" 2 s" 1 at an age of 140 yr. 
For comp arison, the flux rep o rted for Cas A is 3.4 x 10~ 5 ph 
cm" 2 s" 1 (llvudin et al.ll 1994]) . iRenaud et all (120061) examined 
INTEGRAL data from the GC region, acquired within the 
first 2 years of its operations, while searching for 44 Ti-decay 
emission from another claimed young supernova remnant, 
GO. 570-0. 01 8. They showed a map of the Galactic plane 
near the GC with no evidence for emission from Gl.9+0.3. 
They did not detect their target; the limit should apply as well 
to Gl.9+0.3. For an age of 140 yr, their 3a upper limit on 
the yield of 44 Ti is 4 x 10~ 5 M ; at 160 yr, it has increased 
to 6 x 10" 5 M . We believe a deeper INTEGRAL observa- 
tion will provide an important test of the youth of Gl.9+0.3 
and may even help discriminate among progenitor masses (for 
CC) or explosion mechanisms (for la). 

We can make crude predictions of continuum gamma-ray 
fluxes to be expected from Gl.9+0.3, using results quoted in 
Reynolds (2008). The integral photon flux above 100 MeV 
due to the decay of 7r° mesons created by the collisions of 
cosmic-ray ions accelerated in Gl.9+0.3 with ambient ther- 
mal gas is about 6 x 10~ 9 noEsi6 ph cm -2 s -1 , where the explo- 
sion energy is lO 51 ^! erg, and y 0.1 is the fraction of ex- 
plosio n energy in cosmic-ray ions (iDrury. A haroni an. & Volkl 
1994). This should very roughly scale as photon energy E 
or slightly steeper, so that the integrated flux above 0.2 TeV 
is about 3 x 10~ l2 noEsi6 ph cm -2 s -1 , or about 1 mCrab. This 
is about 20ftp times below the H.E.S.S. upper limit in the GC 
region (I Aharonian et al.ll2006h ; with 9 = 0.1 and our density 
estimate of no ~ 10~ 2 cm -3 , we do not expect a detectable sig- 
nal. The very high E m inferred from the synchrotron spectrum 
means that the predicted inverse-Compton spectrum from up- 
scattering cosmic microwave background photons (ICCMB) 
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will be depressed by Klein-Nishina effects. Any TeV upper 
limit can be accommodated by simply raising the magnetic 
field, as it means that fewer electrons are required to explain 
the synchrotron spectrum. But the predictions for ICCMB 
TeV flux are so low that the lower limit on the magnetic field 
corresponding to the current H.E.S.S. limit is less than 1 /iG. 

4. CONCLUSIONS 

We have detected expansion between 1985 and 2007 that 
demonstrates that SNR Gl.9+0.3 has an age of about 100 
years. It also has a synchrotron-dominated X-ray spectrum 
with the highest rolloff energy yet observed, but consistent 
with this very young age. Synchrotron domination requires 
both high shock speeds, for acceleration to energies of order 
100 TeV and above, and low ambient densities, to suppress 
thermal X-ray emission. Both these conditions appear to be 
present for Gl .9+0.3. The implied lack of deceleration is con- 
sistent with the requirement of a lack of thermal emission. 

Our results invite many predictions. Our first, that a new 
radio observation should immediately show expansion, has 
already been verified. Another Chandra observation should 
show it as well. Time variations in the nonthermal fluxes, 
both radio and X-ray, are likely and may give information 
about the density profiles of ejecta and/or ambient medium. 
The INTEGRAL satellite will observe Gl.9+0.3 shortly, test- 
ing the prediction for the 1.16 MeV gamma-ray line. The 
high absorbing column makes optical detection hopeless, but 
the implied high dust density suggests that a search for in- 
frared light echos, such as those recently detected from Cas A 
dKrause et al.ll2005h iKim et al.l 120081: iDwek & Arendtl 12008). 
should be undertaken. At a projected distance of only 300 pc 
from the GC, Gl.9+0.3 is likely to be near substantial local 
dust. A more extensive X-ray observation to search for any 
thermal emission is crucial, as such emission could corrobo- 
rate this picture and also give information on a stage of SNR 
evolution never before observed. 
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